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ABSTRACT 



Aims. Surface brightness asymmetries are a very common feature of stars. Among other effects they cause a difference between the 
projected barycentre and photocentre. The evolution of those surface features makes this difference time-dependent. In some cases, e. 
g. for supergiant stars, the displacement can be a non-negligible fraction of the star radius R, and if R > 1 AU, of the parallax. 
Methods. We investigate the impact of surface brightness asymmetries on both the Gaia astrometric solution and the data processing 
flow with a theoretical approach. 

Results. We show that when the amplitude of the displacement is comparable to the epoch astrometric precision, the resulting astro- 
metric solution of a genuine single star may be, in some cases, of low quality (with some parameters up to lOtr off). In this case, we 
provide an analytical prediction of the impact of the photocentre motion on both^ 2 and the uncertainty in the astrometric parameters. 
Non-single star solutions are found, if allowed for the closest stars. A closer look at the parameters of the orbital solutions reveals 
however that they are spurious (since the semi-major axis is smaller than either its error or the stellar radius). It is thus possible to 
filter out those spurious orbital solutions. Interestingly, for the stocastic solutions, the stochastic noise appears to be a good estimate 
of the photocentric noise. 
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1. Introduction 



The Gaia mission (Perryman et al. 2001; Lindegrenetalj[2008) 
will lead to the determination of the astrometric parameters of 
one billion objects with unprecedented precision. During its five 
years of satellite operations, every detected star will be observed 
on average 80 times, thus yielding its position (an, 5n) at the 
reference epoch, its parallax (m), and its proper motion (/v = 
Ha cos 6, fig). Additional parameters, such as the size of the orbit 
in the case of a binary, can also be adjusted. All these parameters 
derive from a x 1 fit- To allow the comparison between models 
(with diff erent number of paramet ers) the F2 statistics will be 
evaluated (IWilson&Hilfertvlll931h : 



|9v 




(1) 



where v is the number of degrees of freedom (the number of 
measurements minus the number of parameters of the fit). If x 1 
does follow a chi-square distribution with v degrees of freedom, 
then F2 follows a A^(0, l)-distribution irrespective of v. A fit 
with a F2 larger than, say, 3 is considered poor. 
Phenomena such as large convective cells or magnetic spots 
cause time-varying brightness asymmetries on the surface of 
stars. Their presence induces a time-dependent displacement of 
the photocentre with respect to the barycentre t hat possibly alters 
the astrometric parameters derive d by Gaia (Bastian & Hefele 
l2005HEriksson & Lindeprerj2007h . lSvensson & Ludwigl (12005) 
studied the dependence of the photocentric variability on the 
stellar parameters across the Hertzsprung-Russell diagram, from 



white dwarfs to red giants, using CO s BOLD hydrodynamical 
models. They found a linear relation between the standard de- 
viation in the photocentre displacement (o~p) and the pressure 
scale height at optical depth unity. 

There is only one red-s upergiant hydrodynam ical 3D model 
available in literature, dChiavassa et alj [2009), that predicts 
a standard deviation for the photocentre motion driven by 
convection of about 0.1 AU. This value, however, might 
be underesti mated because o f the approximations intrinsic to 
the model. IChiavassa et al] (1201 lh extended the analysis of 
ISvensson & Ludwia (120051) with this red supergiant model. 
They argue that a radical change in the convective pattern is ex- 
pected between the giant and supergiant stages but the large gap 
in T e ff lg between the two available models does not allow a law 
for cr P as a function of T^/g to be unequivocally constrained. 
Hence, they proposed two different laws and performed a first 
analysis of the impact of this photocentric noise on the parallax, 
estimating the number of supergiants and giants for which it may 
be relevant. 

The aim of this work is to investigate and provide an analytical 
prediction of the effects of these time-varying surface brightness 
asymmetries on the astrometric parameters derived by Gaia and 
to study their impact on the Gaia astrometric d ata processing. 
On the basis of the 3D supergiant model from IChiavassa et al.l 
(2009), we simulated the photocentre motion due to convective 
cells and added it to Gaia simulated data (Sect. 2). We discuss the 
results of the astrometric fit of those data. In Sect. 3, we restrict 
our study to the single-star astrometric fit. Only five parameters 
(an, 6q, m, /v, n$) are derived in this solution. In Sect. 4, more 
complex, non single-star solutions are studied. 
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2. The model 

Very large displacements of the photocentre are found in red 
supergiants, which are caused by the presence of extremely 
large convective cells. The flow of hot material towards the sur- 
face and the inflow of cold material around large convective 
cells produce huge surface brightness asymmetries. Evidence of 
this phenomenon are provide d by interferometric observations 
of supergiants shch as a Ori (Buscher et al.ll l990t IWilson et al.l 
19921: iTuthill et al.lll997t TYoung et al.ll200a iTatebe et al.ll2007t 
Haub ois et al.l2009i;IOhnaka et al]201 lb . A numerical 3D model 
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Fig. 1. Snapshot of the 3D s upergiant model in the Gaia G band 
from IChiavassa et aD ((2009). The colour scale denotes different 
temperatures. The 12 M star is modelled with an equidistant 
numerical mesh of 235 3 grid points and a resolution of 8.6 R G 
(or 0.040 AU). The luminosity average over spherical shells and 
over time (i.e., over 5 years) is L = 93 000 ± 1300 L , T e ff = 
3490+13 K,R = 832 + 0.7 R Q , and logg = -0.337 ± 0.001. 



of a red supergian t was developed bv IChiavassa et all (T2009) 
using C0 5 BOLD dFrevtag et al.1 12002: Frevt ag & H6fneil l2008). 
This model was used to predict spectra and intensity maps in the 
Gaia G band for the whole simulation time sequence, namely 
« 5 years with time steps * 23 days apart. A snapshot of the star 
is shown in Fig.Q] 

The projections P x and P y of the photocentre position on either 
axis are central to the present study. We assume that P x and P y 
are random variables (of zero mean) that are both Gaussian and 
MarkoviarQ. Therefore the process correlation function 

C,, v (Af) = E [(P x , y (t) - /-.|/',,|)( /',,(/ + Af) - E[P Xty ])] (2) 

(w here the ope rator E [ ] denotes the mean value over f) reduces 
to (iDoobll 19421) 



C(Af, t) = <T 2 P exp [-Af/r]. 

The random variables P A V at time f, are given by 

P XJ (ti) = exp(-Af,/r)f , A , v (f / _ 1 ) + gi 



(3) 



(4) 



where Af,- = t-, - f,_i, g, is extracted from a Gaussian distribu- 
tion Af(0, s) with s 2 = [1 - exp(-2Af/r)] cr 2 p . Thus, all statistical 
properties of such a Markov chain are determined by the three 

parameters of its mean E [Px,y\ (which we here assume to be 

zero), its variance cr 2 p , and its relaxation time r, i. e. the time af- 
ter which the correlation function becomes "negligibly small", 
which in the present case is cr 2 /e. 
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1 A discrete process is Markovian when its future depends only on its 
present state and not on its past. 
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Fig. 2. Upper panel: photocentre positions extracted from the 
snapshots over the five years of simulation. Lower panel: photo- 
centre positions generated with a Markov chain with cr P = 0.09 
AU and t = 4 x 10 6 s. Points are 23 days apart in both panels. 



We computed the correlation function for the sequence of pho- 
tocentre positions extracted from the snapshots of the 3D model 
in the Gaia G band shown in the top panel of Fig. [2] The re- 
sult is shown as thick lines in Fig. [3] Because of the time step 
of 23 days between two successive snapshots, we are only able 
to evaluate the correlation function on a discrete mesh with a 
23-day time step. We see in Fig. [3] that the correlation function 
falls to zero for At « 60 days. Assuming the correlation func- 
tion to be exponentially decreasing as in Eq. (O, t ^ 46 days 
^ 4 x 10 6 s. The variance cr 2 p is taken as the arithmetic mean 
of cr 2 p t and cr 2 Py , the values of the correlation function at Af = 0: 
crp - 0.09 AU. To check whether the assumed exponential corre- 
lation function is compatible with the one extracted from the 3D 
model data, we generated photocentre positions with cr P = 0.09 
AU and r = 2 - 4 x 10 6 s. An example of the generated posi- 
tions is shown in the lower panel of Fig. [2] We then computed 
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Fig. 3. Thick lines: correlation function of the projections of the 
photocentre positions extracted from the 3D model. Thin lines: 
correlation function of simulated photocentre positions gener- 
ated with an exponentially decreasing correlation function (see 
Eq.0 with crp 0.09 AU and r 2 - 4 x 10 6 s sampled every 
23 days. 
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Fig. 4. Histograms of the P x (left panel) and P y (right panel) po- 
sitions extracted from the snapshots of the 3D model (red striped 
histograms) and the ones generated with the Markov chain (blue 
solid histograms). The actual P x and P y are shown in Fig. [2] The 
black solid line represents the Gaussian distribution used in the 
Markov chain to generate the photocentre positions. 



the correlation function with time-steps of 23 days. The results 
for some series of generated positions are shown as thin lines in 
Fig. [3] and are in good agreement with the correlation function 
for the 3D model points. Furthermore, it shows that we cannot 
distinguish between r = 2 x 10 6 s and r = 4 x 10 6 s. 
we note that the two parameters cr P and r depend on the 
surface gravity and o n the effective temperature of the star 
dChiavassa et al.ll2011l) . This means that the values we derived 
here are only valid for stars with physical parameters similar to 
the ones of the 3D model. In the analysis below, we also consider 
stars with different crp and t in order to provide predictions for 
a wider class of objects. 

In Fig. |U we compare the distributions of P x and P y extracted 
from the snapshots of the 3D model and the ones generated with 
a Markov chain. The solid black line is the Gaussian distribu- 
tion A^(0, <r P ) from which the generated positions are extracted. 
The 80 measurements considered here are not numerous enough 
to clearly reproduce the Gaussian distribution, though it does 
so if we simulate more photocentre positions with the Markov 
chain. We could not verify whether the distribution of the posi- 
tions from the 3D model becomes Gaussian for longer times or 
with more measurements because no additional points are avail- 
able. 

The next step is to add these photocentre displacements gener- 
ated with the Markov chain to simulated Gaia transit data. 
The stellar population considered is a sample of =s 5000 red su- 
pergiants distributed around the Galactic plane (0 < I < 180 and 
-20 < b < 20 in Galactic coordinates). These stars have ab- 
solute V magnitudes between -7 and -3.5 and apparent G mag- 
nitudes between 5.6 and 20. Almost no stars are found at less 
than 1 kpc (for example see Fig . O ■ For such close distances, su- 
pergiants tend to appear brighter than the CCD saturation limit 
at G = 5.6 and will not be obse rved. The sam ple was gener- 
ated with the Besancon model dRobin et alJl2004h an d virtually 
obser ved with the Gaia Object Generator v7.0 (GOG; Isas i et alj 
120101) . We then added to these "observations" photocentric shifts 
generated with a Markov chain. In this paper, we always add to 
all the stars in our sample a photocentric noise generated with 
Markov chains of the same crp and t regardless of the physical 



parameters of the star. 

The conclusions we draw in the analysis below shiuld be inter- 
preted in a statistical sense, as star-by-star effects are closely re- 
lated to the particular Gaia scanning law and actual photocentre 
configuration. 

3. Impact on the single-star parameters 

As mentioned in the previous section, the goal of our work is to 
study the impact of surface brightness asymmetries on the stel- 
lar astrometric parameters derived by Gaia. The observations, as 
described in the previous section, are obtained by adding pho- 
tocentre motions generated with a Markov chain to Gaia transit 
data. They were fitted with the five-parameter single-star model 
(qtq, 60, m, fi a -, n$). In Fig. [5] we show some results concerning 
the parallax. The upper panel of Fig. [5] illustrates the difference 
Am between the parallaxes determined with (m a ) and without 
(ifrna) the surface brightness asymmetries, in units of the statisti- 
cal error in the parallax cr w . We note that cr m will be provided as 
output in the Gaia catalogue and is derived from the covariance 
matrix of the parameters obtained in the data fitt ing process. As 
shown in Appendix A of IChiavassa et al.l d201 ll) and of this pa- 
per, the quantity cr w is independent of the presence of surface 
brightness asymmetries. The effect of the photocentre displace- 
ment can be quite significant as it can induce deviations in the 
derived parallaxes of up to 10cr for d ~ 1 kpc. Nevertheless, we 
not that all effects become negligible for stars farther out than 
5-6 kpc. At this distance, the epoch astrometric precision, which 
is about 0.03 mas for a star brighter than G = 12, is comparable 
to the angular dimension of the added photocentre displacement 
(0.09 AU on average). A smaller photocentre noise will clearly 
have an impact only for more nearby stars, and the opposite for 
larger crp. For example, at cr P - 0.02 AU, only the most nearby 
stars (d < 2 kpc) have a Am greater than 3cr m . On the other hand, 
at crp = 0.2 AU, a deviation up to 30<r is reached. 
The expected value of Am can be easily computed (details in 
Appendix A) assuming that the photocentre displacements are 
taken from a Gaussian distribution with variance crp. We note 
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Fig. 5. Upper panel: relative variation in the single-star parallax 
{(fijna ~ &a)/o- m ) resulting from the addition of surface bright- 
ness asymmetries. Lower panel: same as upper panel but with 
the standard deviation predicted from Eq. (|5]l added (red points). 
The parameters of the Markov chain used to generate the data 
are cr P = 0.09 AU and t = 4 x 10 6 s. 



that with this approach the positions of the photocentre are 
treated as uncorrected, which, as we have seen, is not strictly 
true. The mean value of Am is zero because the photocentre posi- 
tions, on average, coincide with the geometric centre of the star. 
This is also true for correlated photocentre positions if enough of 
them are considered. On the other hand, cr^, the standard devi- 
ation in Am, indicates how much individual m a and m na values 
can differ: 



N- 



<T„ + CT 



(T- 



(5) 



where d is the distance in kpc, N is the number of measurements, 
cr w is the epoch astrometric precision , <Xp and cr^ are the vari- 
ances in the parallax factors dGreenl ll985. and Eqs. IA.191 and 
IA.201 i. In the lower panel of Fig. [5] we have plotted the cr^ val- 
ues estimated from Eq. (O on top of the Am resulting from the 
single-star astrometric fit. The agreement between the two sets 
of points is remarkably good, meaning that cr^m can be used as 
a realistic estimate of the impact of the surface brightness asym- 
metries on the accuracy of m a . 



The Gaia end-of-mission astrom etric accuracy cr^ for a 
feet" single star can be written as JdeBruiine 200ll) 



per- 



(6) 



In this expression, m is an overall, end-of-mission contingency 
margin (m = 1.2 for a 20% margin), g w is a dimensionless geo- 
metrical factor relating to the scanning geometry, cr w is the single 
measurement error (including all the error sources), and N is the 
number of transits. We can rewrite Eq. (0 for cr A ^ in terms of 
Eq. © 



m 2 gj (Tp 



Pa 



' Ps cr© 
oc 



(7) 



We see that cr^w is only related to the angular size of the photo- 
centre displacement induced by the surface brightness asymme- 
tries (erg = 2p) and the scanning geometry. There is no depen- 
dence on the instrumental sensitivity cr K because cr AS . is intrinsic 
to the physics of the star. On the other hand, the capability 
of Gaia to detect the noise induced by the surface brightness 
asymmetries, does depend on cr w such that 



0"C7 0\.< 



(8) 



For red supergiants, cr w is roughly constant because G < 12.6 in 
most of the cases. 

Finally, if we consider ^ M -, which is the parallax inaccuracy (ow- 
ing to the surface brightness asymmetries) relative to the paral- 
lax, we find that it is independent of the distance of the star 



OC (Tp. 



(9) 



We can now compute 2** and 2^s. for a red supergiant star us- 
ing mean values for the scanning- law parameters _ 1.91, 
m = 1.2, N = 78, cr M , = 30 //as dde Bruiindl2005t iLindegrenl 
120101) . and crp = 0.09 AU. For supergiants located at 1 kpc 
(ere = 0.09 mas), Eq. © predicts that < 4.3 for 68% of 
the sample and values larger than 13 for 0.3% of the sample. For 
supergiants at 5 kpc (cr = 0.018 mas), these numbers are 0.8 
and 2.4, respectively. Similarly, Eq. <j9j predicts that < 3.4% 
for 68% of the supergiants (independently of their distance) and 
values larger than 10.3% for 0.3% of the sample. 
A similar behaviour is found for the four other parameters (al, 
Sq, yL a -, and Expressions equivalent to Eq. © are derived 
in Appendix A. we note that for position and proper motion, 
we find values that are up to 15<x off, while the photocentric 
noise was generated with the same Markov chain parameters as 
in Fig. 

We now analyse the accuracy m, rue - m a (as opposed to Am = 
m a - m„ a ) of the "fitted single-star" parallax with respect to the 
true parallax. In the three panels of Fig. [6] we plot this differ- 
ence, normalised by the estimated error in the parallax cr m , for 
different slices in distance. For stars far away (right panel), the 
distribution is clearly A^(0, 1) as expected. As we consider more 
nearby stars, the distribution becomes wider and wider. For stars 
in the range 3.5 < d < 5.5 kpc (central panel), the histogram 
is approximated well by a A^(0, 2.5) distribution, and for stars 
with d < 3.5 kpc (left panel) by A^(0,4). If we increase crp, 
we obtain wider and wider kpc distributions. For cr P = 0.2 AU, 
only stars between 7.5 kpc and 8 kpc follow a A^(0, 1) distribu- 
tion. On the other hand, lowering crp produces narrower distri- 
butions. For crp = 0.02 AU, the distribution becomes A^(0, 1) for 
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Fig. 6. Histogram of the relative error in the single-star parallax 
((j^mie - & aliens) resulting from the addition of surface bright- 
ness asymmetries for supergiant stars at different distances. The 
distribution for stars in the right panel (d > 5.5 kpc) closely 
reproduces the expected A^(0, 1). In the central panel (3.5 < 
d < 5.5 kpc), the distribution is wider with a best-fit N(Q,2.5). 
Finally, for the left panel (d < 3.5 kpc) the cr of the distribution 
increases to four. Parameters of the Markov chain used here are 
a-p = 0.09 AU and t 
are at the same scale. 



4 X 10 s. The x-axes in the three panels 



stars at all distances indicating that at this level of photocentric 
noise, no effects are expected in Gaia data processing. We recall 
that actual stars have cr>, which is dependent upon their abso- 
lute magnitude, but th e exact relationship is not known so far 
dChiavassa et al.ll201ll) . 

As already mentioned, we have used a five-parameter-model 
for our set of supergiants because we know we are dealing with 
genuine single stars. However, the F2 values for these solutions 
(Fig. [7] top panel) can exceed the threshold value three by more 
than one order of magnitude. An estimate of the increase in F2 
induced by the presence of surface brightness asymmetries can 
be easily obtained by noticing that Xa> which is the x 1 m the 
presence of surf ace brightness asymmetries, can be written as 
dChiavassa et al.ll201 ll Eqs. 11 and 12, and Appendix B of this 
paper) 



' X n 



+ N 



a-p 



(10) 



where xl a is the x 1 of the fit of the data without surface bright- 
ness asymmetries, N is the number of measurements, d is the 
distance in kpc, and cr w is the Gaia single-measurement error. 
This expression can be easily converted into the analogous ex- 
pression for F2 



F2 a — F2 na + J ^ 



1/2 



cr 2 V /3 
1 + ^1 

U W } 



(11) 



This quantity is plotted in the lower panel of Fig. [7] on top of 
the results obtained for the single-star solution. There is on aver- 
age remarkable agreement. As we previously pointed out, a case 
by case agreement is not to be expected. When we change the 
parameters of the Markov chain and set o> = 0.2 AU, F2 a in- 
creases to 100 for more nearby stars. No stars at d < 2.5 kpc 
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Fig. 7. Upper panel: F2 for the single-star fit of the simulated 
data with the addition of surface brightness asymmetries. A 
F2 > 3 indicates a non-satisfactory fit. Lower panel: the same 
as the upper panel but with the theoretical prediction (Eq. ( TPTl i) 
added on top (see text). 

have F2 a < 3, and we find values larger than three at all dis- 
tances. In contrast if we look at what happens with cr P = 0.02 
AU, F2 a becomes larger than three only for a negligible fraction 
of the stars, and it never exceeds seven. 



4. Non-single-star processing 

We have so far treated the supergiants in our sample as single 
stars, as we know that they actually are. Nevertheless, when pro- 
cessing real Gaia data, we will not know in advance the kind 
of object we will be analysing. When we do not restrict the 
data processing flow to the five-parameter solution, stars with a 
single-star F2 larger than three are further processed. A chain of 
solutions of increasing complexity is applied until F2 becomes 
smaller than three. In the acceleration solution, a possible vari- 
ation of the proper motion is added to better fit the data. In the 
orbital solution, the parameters of the orbit of a binary system are 
added. If neither the single-star solution nor the non-single-star 
solutions (acceleration or orbital) are acceptable on the basis of 
their F2, a stochastic solution is computed: some noise is added 
to the along-scan measurements in order to lower to zero the F2 
of the single-star solution. 

For our sample, an acceleration solution, an orbital solution, or a 
stochastic solution are found for some fraction of the stars, and 
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this fraction depends on the parameters of the Markov chain. 
The distribution of solutions among these possible classes is 
plotted in Fig. [8] for three different sets of crp and r. The num- 
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Fig. 8. Distribution of the different types of solutions found as a 
function of the distance. In the three panels are shown the results 
for the same sample of stars to which different photocentric mo- 
tions have been applied, as indicated by the corresponding crp 
and t in the figure label. 



ber of single-star solutions decreases drastically as we increase 
crp. Moreover, the impact on the astrometric parameters then in- 
creases, implying that the number of solutions in which they are 
unaltered decreases. We see as well that all the nearby (d < 4 
kpc) single-star solutions disappear. At the same time, we see an 
order-of-magnitude increase in the number of stochastic solu- 
tions. Finally, we point out the decrease of the fraction of orbital 
solutions as we decrease r. This is a consequence of the loss of 
correlation between temporally close measurements that makes 
it harder to fit the points with an orbit. 

That single-star solutions are rejected by the Gaia pipeline even 
though we are dealing with genuine single stars is unfortunate. 
Hence, we analysed these solutions to see whether it is possible 
to identify them and filter them out. We focused our attention on 
the orbital solutions as they can represent a significant fraction 
of the solutions for our sample. 

As a first step, we analysed the significance of the orbital so- 
lution by looking at the ratio a/cr a , where a is the semi-major 
axis of the orbit. If a is smaller or comparable to its error cr a , the 
solution is not significant as this is an indication of null orbital 
parameters. In the left panel of Fig. [9] we show the ratio a/cr a 
for a Markov chain with cr P = 0.2 AU and t = 4 x 10 6 s and 
conclude that a large fraction of these solutions has a/cr a < 1 . 
They are thus spurious and may easily be flagged as such. For 
these solutions, the reduction in F2 under the acceptance value 
of 3 is caused by the x 1 decreasing faster than the number of 
parameters increases. As a second step, we studied the stars for 
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Fig. 9. Left panel: distribution of the ratio a/cr a , with a the semi- 
major axis of the orbit and cr a its uncertainty. Solutions with 
a/cr a < 2 are not significant and may be rejected. Right panel: 
ratio a/R, with R the stellar radius for stars with a/cr a > 2. All 
orbits appear to have a semi-major axis smaller than the stellar 
radius and can thus be rejected. Markov chain parameters: crp = 
0.2 AU and t = 4 x 10 6 s 



which the orbital solution is significant (i.e. located in the peak 
at about a/o~ a = 6 in the left panel of Fig.|9]l and we computed 
the ratio a/R, where R is the radius of the star. This quantity is 
plotted in the right panel of Fig. [9] and is always smaller than 
unity. We may thus reject all those solutions as they are unphys- 
ical, representing orbits with separations smaller than the radius 
of the stars under study. Moreover, we note that the size in AU of 
the semi-major axis found (Fig. ITOl is always of the same order 
of magnitude as the photocentre motion added with the Markov 
chain, 0.2 AU being about the minimum value found. 
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Fig. 10. Semi-major axes a of the orbits derived with the orbital 
solution. The physical size of a in AU is comparable to the pho- 
tocentre motion added to the simulated Gaia data. Markov chain 
parameters :o> = 0.2 AU and t = 4 x 10 6 s 



Finally, we discuss the stochastic solutions. Stochastic noise 
is meant to represent any source of "noise" found in the data that 
cannot be modelled with any of the other models. We expect it 
to be correlated with the noise we added to the single-star Gaia 
data. A histogram of the stochastic noise (in AU) is shown in 
Fig.QT|for three different sets of parameters of the Markov chain. 
We show examples with extreme parameters in order to more 
clearly illustrate the occurring effects. In the left panel, we show 
the results for simulated photocentric motions with cr P = 0.2 
AU and t = 4 x 10 6 s. The blue horizontally striped histogram 
represents the distribution of the stochastic noise obtained by 
the data processing. It is clearly neither symmetric nor centred 
at the value of cr P - 0.2 AU as expected (vertical dashed line). 
This is because in the fitting process, the position of the star is 
assumed to be at the mean photocentre position, which does not 
necessarily coincide with the geometrical centre of the star (ex- 
ample in the top panel of Fig. [2] where the mean photocentre po- 
sition is in the first quadrant rather than at the origin). Once we 
quadratically add to the stochastic noise the difference between 
the mean photocentre position and the true position of the star, 
the resulting distribution (green diagonally striped histogram) 
becomes centred at 0.2 AU as expected. This positional shift and 
the asymmetry of the distribution is caused by the length of the 
Gaia mission (5 years) not being long enough for the distribution 
of photocentre positions to be Gaussian. Moreover, that mea- 
surements are taken at time intervals shorter than the relaxation 
time of about 40 days and the non-uniformity of the time sam- 
pling play a crucial role. The minimum time between two con- 
secutive Gaia observations of the same target is about two hours. 
When we decrease r to 4 x 10 3 s (i.e. about one hour, meaning 
that all measured photocentre positions are almost uncorrected), 
the histogram becomes almost symmetric and centred on 0.2 AU 
(central panel of Fig. fTTb - If we keep t = 4 x 10 6 s and decrease 
crp to 0.08 AU (right panel of Fig.lTTb. the stochastic noise distri- 
bution is still asymmetric and not peaked at 0.08 AU. However, 
the effects are mainly smaller than for the 0.2 AU case, merely 
because of the smaller size of the photocentric shifts generated. 

From this discussion, we can conclude that the stochastic 
noise can be used as an estimate of the photocentric noise present 
on the surface of the star. This correspondence is at its closest for 
short correlation times and small photocentric motions. Under 
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Fig. 11. The blue horizontally striped histograms show the 
"stochastic" noise added to Gaia measurements in the compu- 
tation of a stochastic solution. The green diagonally striped ones 
represent the stochastic noise plus the offset of the mean photo- 
centre position with respect to the centre of the star. The three 
panels show results for three different sets of parameters of the 
Markov chain. The vertical dashed line corresponds to the stan- 
dard deviation crp in the photocentre displacement. 

these circumstances, the distribution will be more symmetric and 
will peak at the correct value. 



5. Conclusions 

Surface brightness asymmetries are a common characteristic 
of stars and cause time-varying positional shifts between 
the photocentre and the barycentre. We have analysed their 
effects on the parameters derived from Gaia simulated data for 
supergiants. To achieve this goal, on the basis of a 3D model of 
a supergiant, we have simulated photocentric displacements to 
be added to along-scan measurements. We have then performed 
a five-parameter fit. For such a single-star solution, we have 
shown that the deduced parallaxes can be up to I3cr m off for 
a supergiant at 1 kpc with cr P = 0.09 AU and 2M. U p to 10%. 
We have shown as well that 2M decreases with distance as 
-p-, but in contrast, 2M does not vary with distance, being 
proportional to crp. At the same time, the quality of the fit can 
be very poor, with F2 parameters exceeding the threshold value 
of 3 by more than one order of magnitude. If more complex 
solutions are searched for, even though F2 becomes smaller 
than three, there is no improvement in the accuracy of the 
derived parameters and, at the same time, the computational 
effort is substantially increased. Spurious binary solutions can 
be rejected on the basis of significance tests and the comparison 
of the size of the semi-major axis with the radius of the star. 
When a stochastic solution is found, we observed a correlation 
between the stochastic noise and the Gaussian noise injected 
in the Markov chain to simulate the photocentric motion. This 
information could be exploited to obtain, at least, the order of 
magnitude of the convective noise present. 
From the standpoint of the Gaia data processing, if a star could 
be identified as a supergiant independently of its parallax, 
a 'supergiant flag' may then be used as a warning that the 
astrometric parameters may be far off. Classification methods 
independent of the parallax include the spectroscopic determi- 
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nation of the atmospheric parameters (T e ff, logg) from the Gaia 
spectrometer, or the classification as an irregular variable from 
the repeated photometric measurements. Irregular variables 
exist in two flavours, ho wever: supergiants a nd stars at the base 
of the red giant branch dJorissen et alj|1997h . The two classes 
of irregular variables differ by several magnitudes in their 
respective absolute magnitudes, so that a 10% relative error in 
the parallax of the supergiant caused by the surface brightness 
asymmetries is not large enough to intermingle the two classes 
and introduce confusion. 

To identify true binaries among supergiants despite the con- 
fusion introduced by the surface brightness asymmetries, an 
indirect way is provided by the detection of composite colour 
indices (in addition to the orbital separation probably now being 
significantly larger than the stellar radius). Such composite 
colour in dices for systems i nvolving supergiants were com- 
puted by Mal kov et al.l d201 ll) . with the conclusion that these 
composite indices differ by at least 0.3 mag from the nearest 
single-star indices in the best two-colour diagram, thus allowing 
the identification of the binary systems. 

This complementary non-astrometric information should thus 
allow for a more efficient and reliable evaluation of the sig- 
nificance of the astrometric solution in the presence of surface 
brightness asymmetries. 
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Appendix A: Estimation of the accuracy on the 
astrometric parameters 

Gaia will provide N measurements (w) of the along-scan posi- 
tion of each star and hence the astrometric solution will be com- 
puted by solving an overdetermined system of N equations for 
/?„„, the n unknowns (a* () , So, zo~, n a ",Hs in the case of a single-star 
solution discussed here): 



Xp na + e = w, 



(A.l) 



where e is the vector of errors that we assume to be random 
and mutually uncorrected, i.e. they have zero mean and their 
covariance matrix W is diagonal. We define X to be the matrix 
of known coefficients 



X 



sin f?i cos 0i f Wl ?isin#i ficos^i 
sin 02 cos 82 f W2 ?2 sin 62 ?2 cos 62 



sin 6m cos 9^ fw N ?w sin 9^ cos 9, 



(A.2) 



,v ) 



where 9 is the position angle of the scan, t is the observing time, 
and f w = r(P a sin 9 + P$ cos 9) is a lin ear combination of the par- 
allax factors P a and P$ dGreenll 19851) . and r is the radius vector 
of Gaia's orbit at time t in AU 

P a — cos € cos (x sin © — sin c? cos O, (A. 3) 

Pg = (sin e cos 6 - cos e sin a sin 6) sin O - cos a sin S cos 0, 

where O and e are, respectively, the longitude of the Sun and the 
obliquity of the ecliptic, both at time t. 

We wvaluate /}„„, the best estimate of the parameters /3 na in 
Eq. dA.U . by minimizing;^ 2 



X 2 na = (w- xp na ) T W \w - XJ3 na ), 
which can be expressed by 



frna = {X T W- 1 X)-\X T W- 1 W). 



(A.4) 



(A.5) 



The difference between the true parameters and the estimated 
ones is easily obtained using Eqs. ( IA.5t and (IA. lb 

Pna -Pna = (X T W~ X XT' X T W~ X €, (A.6) 

as well as Z, the covariance matrix of the parameters 

Z = E[0 na - (J3 na ))0 na - (J3 na )) T ] = (X T W-'XT 1 . (A.7) 

We note that Z is diagonal because all the off-diagonal terms 
contain sin 9, cos 9, or t factors that averages to zero. The mean 
parallax factors are also zero because of their dependence on O. 
We now consider a star with time-varying surface brightness 
asymmetries. Their presence causes a shift in the photocentre 
position and a consequent shift (Aw) in the observed along-scan 



E. Pasquato et al.: Limits in astrometric accuracy induced by surface brightness asymmetries in red supergiant stars 



9 



position with respect to a star with a uniform disk. 
This can be expressed as 

Aw/ = — [a; cos 9j + bi sin 6{\ = —S>u 
a a 



(A.8) 



where (a,-, bi) are the coordinates of the photocentre in a refer- 
ence system with its origin at the centre of the star (expressed in 
AU) and d is the distance of the star in kpc. We assume that a, 
and bi follow a Gaussian distribution with zero mean and stan- 
dard deviation cr P . 

The astrometric problem to be solved now to get the parameters 
p a is analogous to Eq. ( IA. lb 



Xf) a + € - W + Aw. 



(A.9) 



The estimate of the parameters is again found by x minimiza- 
tion 



xl = W- 1 [(w + Aw) - Xp a ] T [(w + Aw) - Xpa\ 
and the result, in analogy with Eq.( lA.51 >. is 
j} a = (X T W- 1 Xy 1 (X T W- 1 (w + Aw)). 



(A. 10) 



(A.ll) 



With the use of Eqs. (IA.5t and dA.l U , we can produce an ex- 
pression for the effect of the surface brightness asymmetries on 
the parameter estimates 



Pa -Pna = (X' W^XyHX 1 W~'Aw) = Z(X' W~ l Aw). (A.12) 

For a non-variable star, the elements of Z and cr w ., will be ap- 
proximately all the same as they are a function of the apparent 
magnitude G and we can rewrite Eqs. ( lA.5b . ( IA.7I ). and ( IA. 12b 
for the derived parameters and the covariance matrix, as 



Pna = (X T X)-\X T W), 

p a = (X T Xy\X T {w + Aw)), 
Z = o-l{X T X)-\ 



(A.13) 
(A. 14) 
(A. 15) 



We note that the covariance matrix Z is independent of the pres- 
ence of the brightness asymmetries, which means that the un- 
certainty in the deduced parameters will be the same. That is 
because what enters the covariance matrix are basically the mea- 
surement errors. 

In reality, the accuracy (to be distinguished from the precision) 
of the parameters can be quite different in the two cases. Indeed, 
the parameters derived in presence of brightness asymmetries 
will be different from the ones found for a uniform-disk star, 
their difference being given by: 



hfi - Pa~ Pna = Z(X T W~ l \w) = ~ 

o-ld 



value of the derived parameters caused by the presence of the 
brightness asymmetries. For the parameter Pj, this means that we 
need to estimate the variance in N{XijDi).. We can easily com- 
pute those factors for all the parameters where (a) = (b) = 0, 
o-\ = cr\ = a% and <f> = 0. 

For the first two parameters al and 6, we have to compute the 
variance of {XjjDi) i = (sinff)D and {XjjD^ = (cos 8)D, which 
leads to 



and 



N °% ay 



2crl d' 



(A. 17) 



(A.18) 



For the parallax (XyD,). = f v D, we need the variance in the 
parallax factors, which is easily obtained by noting that the only 
term that is not constant in time is O and neglecting the variations 
in r. 

cr^ = -(cos 2 ecos 2 a + sin 2 ff), (A. 19) 

cr 2 4 = - ^(sin e cos 6 - cos e sin a sin S) 2 + cos 2 a sin 2 tf] . (A. 20) 



Hence 



0~ Am — O* ib a -m m 



+ a% 



Ncr^ay 
2 o-l d ■ 



(A.21) 



Finally for the two components of the proper motion fi a * and 
ps, we have to compute the variance in {XijDi) j = f(sin#)D and 
(XijDi). = f(cos ff)D, which leads to 



0"A/j„. - ~pL a >- a -pL a >- m - 



and 



2 IN °~ n„, (Tp 



' V 2 o-l d 



, 2 Na^o-p 



(A.22) 



(A.23) 



,(A.16) 



where cri are the diagonal elements of the covariance matrix Z, 
X is the matrix of the known coefficients described in Eq. ( 1A.21 i. 
and D describes the photocentre motion (see Eq. IA.8b . In this 
equation, all the terms are known apart from Ylj=\ XjjDj. These 
terms are on average zero because the factors are all uncorrected 
and (Di) = 0. As we expect to see some effects on the deduced 
parameters because of the presence of surface brightness asym- 1 cr 2 

metries, it is interesting to compute cr 2 ., the standard deviation — x\ a + N — I — - ] 

of p a - p na . This quantity gives an estimate of the errors in the 



Appendix B: Expression for x 1 in the presence of 
surface brightness asymmetries 

We estimate the increase in the x 1 caused by the presence of 
surface brightness asymmetries. We can rewrite Eq. (IA. 10b for 
xl in a slightly different form 

Xl = W l [w-Xp na + Aw-XAP)] T [W-Xp na + Aw-XAP)](B.1) 

where Ap is p a - p„ a . The minimum xl will t> e > by defini- 
tion, the one for the deduced parameters, hence with Ap = A/3 
(Eq. ( IA.16b ) and with p„ a = P„ a (Eq. dA.5b ). Using those values, 
we obtain 

xl = W' 1 \(w - Xp na f (w - Xp na ) + Aw T Aw + (XApf (XAp) + 



-2Aw T (XAp) + 2 (w - Xp na f Aw - 2 (w - Xp na f (XAp) 



■(B.2) 



The last four terms are negligible with respect to the first two, 
leading to 



(B.3) 



